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Abstract 

The main purpose of this article is to disseminate among a wide au- 
dience of physicists a known result, which is available since a couple of 
years to the cognoscenti of differential forms on manifolds; namely, that 
charge conservation implies the inhomogeneous Maxwell equations. This 
is the reciprocal statement of one which is very well known among physi- 
cists: charge conservation, written in the form of a continuity equation, 
follows as a consequence of Maxwell equations. We discuss the conditions 
under which charge conservation implies Maxwell equations. The key role 
played by the constitutive equations is hereby stressed. The discussion is 
based on Helmholtz theorem, according to which a vector field is deter- 
mined by its divergence and its curl. Green's functions are also shown to 
play a fundamental role. We present all results in three-vector, as well 
as in tensorial notation. We employ only those mathematical tools most 
physicists are familiar with. 
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1 Introduction 



Maxwell equations are frequently introduced - using gaussian units - in the 
following form pQ: 



V D = 4np (1) 
T ldD 4tt 

VxH -c^ = T j ® 

V B = (3) 
1 <9B 

VxE+ = o. (4) 

C OT 

Equations and © are called inhomogeneous - or Maxwell equations with 
sources -, while @ and (@J are called homogeneous, or source-free equations. 
The four equations constitute a closed system because the couples (D,H) and 
(E, B) are related to each other through the so-called "constitutive equations". 
It is however not unusual to stress the fact that the equations with sources are, 
to some extent, conceptually different from the source-free equations. Indeed, 
Eqs.Q and gj) can be understood as expressing a purely mathematical state- 
ment. To see this we start by considering a scalar field ip(t, x) and a vector field 
A(£,x), which are continuously differentiable but otherwise totally arbitrary. 
Then we construct the fields Be VxA and E = — V<y9- d t A/c. Eq.© is then 
identically satisfied because the divergence of a curl vanishes. If we now take 
the curl of E and use the fact that the curl of a gradient vanishes, we see that 
Eq.|@J also holds true identically. We conclude that Eqs.Q an d Q are satisfied 
by arbitrary fields, as long as these fields are constructed as above, starting from 
the given fields tp and A. These two equations are therefore not characteristic 
of the electromagnetic field. They can be understood as a mathematical state- 
ment telling us that there are fields, ip and A, out of which we can construct 
E and B. The electromagnetic nature of these fields depends on the fact that 
they have to satisfy also equations (JJJ and (0), as long as D = D(E,B) and 
H = H(E, B). Equations and (0) are the ones possessing a truly physical 
content. They are the ones which contain the sources that produce the field. 
It is the particular way in which these sources are related to the fields, what 
makes up the physical content of these equations. 

Let us now turn to Maxwell equations as they are often written in tensorial 
form: 



An 

d,F^ = —f, (5) 

dpFvx + dvFxp + dxFpy = 0. (6) 

Here again, starting from an arbitrary four- vector A^(x) we may define an 
antisymmetric tensor F M „ = d^A u — d v A^. It is easy to see that this tensor 
identically satisfies the homogeneous equation © , which is a Bianchi type iden- 
tity. As before, if our has to describe an electromagnetic field, then it has 
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to satisfy the inhomogeneous equation d^d 11 A v — d v d^A^ = 4irj" /c, which is 
another form of Eq. (JSJ . 

Summarizing, we can say that the homogeneous Maxwell equations can be 
considered as entailing a mathematical statement about the nature of the fields 
E and B, or -correspondingly - about the tensor Fuv- The inhomogeneous 
Maxwell equations in turn are the ones possessing physical content. We must 
postulate that the electromagnetic field has to satisfy them. Now, all these 
things are very well known. What seems to be not so very well known is the 
fact that the inhomogeneous equations by themselves are also not characteristic 
of the electromagnetic field. Indeed, suppose we are given a scalar function 
p(t, r) and a vector function r), both of which go to zero sufficiently rapidly 
as r — > oo, and being such that they satisfy the equation 

d t p + V-j = 0. (7) 

Then there exist vector fields D(t, r) and H(i, r) satisfying the inhomoge- 
neous Maxwell equations 



V D = p (8) 
VxH = j + ftD. (9) 

The existence of D and H can be proved by explicit construction. Such a 
construction rests on Helmholtz theorem [3 El EI 00) which is discussed below. 
For now, it suffices to say that - loosely speaking - "a vector field is determined 
by its divergence and its curl". Thus, according to Helmholtz theorem, Eq.JHJ 
can be solved for D (though the solution is not unique) when p is given. From 
equations © and J7J we see that V-((9 t D + j) =0. Applying Helmholtz theorem 
again we can show that there is a field H whose curl is 9{D + j. This is equation 
JJjJ. Note that we have written the inhomogeneous Maxwell equations in MKS 
units, which are the convenient units for what follows. 

We see then that the continuity equation (JJJ entails the inhomogeneous 
Maxwell equations. The continuity equation expresses the conservation of some- 
thing. This something must not necessarily be electric charge. It could be mass 
as well, or any other quantity - like probability, for instance. We are thus led to 
conclude that the inhomogeneous Maxwell equations are also not characteristic 
of electromagnetism. They hold true whenever something is conserved. Putting 
things this way we bring to the fore the fundamental role played by the consti- 
tutive equations, D = D(E, B) and H = H(E, B), whatever their precise form 
might be. They constitute the link between the homogeneous and the inhomo- 
geneous Maxwell equations. It is this link what turns the four equations into a 
closed system. Neither the inhomogeneous nor the homogeneous equations by 
themselves are characteristic of electromagnetism. They must be linked to one 
another in order to conform a closed system of equations that is characteristic 
of electromagnetic phenomena. 

In the following section we discuss Helmholtz theorem. Although this theo- 
rem can be found in several textbooks and articles, for our purposes it is useful 
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to present it in a form which brings to the fore its connection with Green's 
functions. 



2 Helmholtz theorem 

Here we discuss Helmholtz theorem by following an approach which is slightly 
different from the one presented in several textbooks. Helmholtz theorem states 
that a vector field v is completely determined by giving its divergence and its 
curl, together with its normal component, n • v, at the boundary of the domain 
where such a vector field is to be determined. For physical applications it is 
natural to take as "boundary" an infinitely distant surface and v vanishing there. 
Helmholtz theorem then says that we can write v in terms of two potentials, U 
and C, in the form 

v(x) = VJ7(x) + VxC(x), (10) 

where U and C can be expressed in terms of the divergence and the curl of 
v(r), respectively. Now, put in this form, Helmholtz's theorem might appear 
as a result that is rather awkward to prove. Let us thus try to lay bare what 
motivates it. To this end, consider the following two vector identities, in which 
the Laplacian V 2 appears: 

V-(VU) = V 2 C/ (11) 
Vx(VxC) = V(V-C)-V 2 C. (12) 

Add to these relations the equation satisfied by a Green function G(x, y), 
on which we impose the condition that it vanishes at infinity: 

V 2 G(x,y) = <i 3 (x-y), (13) 

G( "> " -SFjT^r (14 » 

By means of G(x, y) we can introduce U and C as "potentials" associated 
with two given "densities", p and j, through 



We assume that p and j vanish at infinity. The potentials then satisfy 

V 2 L/(x) = p(x), (17) 

V 2 C(x) = -j(x). (18) 

V • j(x) = => V • C(x) = 0. (19) 
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The validity of Eas. l|17ll8|l follows directly from the definitions given by 
Eas. (|15ll6|l . together with Eas. l|13ll4|l . In order to see that V • j(x) = implies 
that C(x) is divergenless we need a little more ellaborated calculation: 



V • C(x) 



1 

4tt 



j(y) 
|x-y| 



^/j(y)-v 5 



47T 

1 

47T 



d 3 y = 
1 



|x-y 



l 

-in 

d 3 y 



j(y) • v x 



|x-y 



d 



3 

i=l y 



1 



-1 

47T 



j(y) • n (y) 
|x-y| 



dS 



3 

£ 

i=l 



0*(y) 



x - y| <9y l 



d 3 ^20) 
(21) 

(22) 
0. (23) 



In the last step - which resulted from an integration by parts - the volume 
integral was replaced by a surface integral using the divergence - or Stokes - 
theorem. Such a surface integral vanishes when the volume of integration goes 
to infinity, because j has been assumed to vanish at infinity. The second term 
vanishes because of the requirement V • j = 0. 

From Eqs.JTTJl and (O together with (|17I18I19|) . we see that 



V-(Vtf) = p(x), (24) 
V x (V x C) = j(x). (25) 

This suggests us to define a field v = VU + V x C. Such a field satisfies 

V-v = p(x), Vx v=j(x). (26) 

This way we arrive naturally at the following statement: if we are given 
the divergence p(x) and the curl j(x) of a vector field v(x) which vanishes at 
infinity, then we can write this field asv = V£/ + VxC, where U and C are 
given in terms of p and j by Eas. (|15|l and (|16|) . In order to see that p and j 
uniquely determine v, it suffices to show that when both the divergence and the 
curl of a field vanish, then the field itself vanishes identically. This follows from 
what we have done so far. Indeed, we have shown that the following equation 
holds true identically: 



v(x) = -V x ^ G(x, y)V • v(y)d 3 y) +V X X (J G(x, y) V x v(y)rf 3 ^ , (27) 

with the Green's function G(x, y) satisfying Eas. (|13I14|I . Hence, if V • v = 
and V x v = 0, then v = 0. We conclude that given two fields, vi and V2, 
having the same divergence and curl, they must in fact be the same field. This, 
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because their difference v = vi — V2 vanishes identically, as a consequence of 
V • v = and V x v = 0. 

Finally, let us first note that Eg. l|27|) holds for Green's functions other than 
the one defined in Ea. l|14|) . Indeed, the only property we need to assume about 
the Green function G(x, y) is that it be of the form G(x — y). This is true 
anyway, whenever G(x, y) fulfills Ea. (|13[l . As to the field v, it has been assumed 
to vanish at infinity. In fact, it suffices to assume that it vanishes faster than 
1/r for large r. Note also that if we prescribe only the divergence V • v = p(x) 
of a field, then what we can deduce from this sole condition is that 



v(x) = -V x (J G(x, y)V • v(y)d 3 j/) + V x x Z(x), 



(28) 



with Z(x) arbitrary. If we instead prescribe only the curl V x v = j(x) of a 
field, then we have 



(29) 



v(x) = V x x G(x,y)V x v(y)d 3 ^ + V x V(x), 
with V(pc) arbitrary. 

3 Maxwell equations and Helmoltz theorem 

We have discussed Helmholtz theorem in the framework of i? 3 . That is, the 
vector fields we have considered are of the form v(x). However, all the results 
we have obtained so far remain valid if we assume these fields to depend on 
a set of additional parameters. They can be assumed to have been there all 
the way, but without having been shown explicitly. Let us denote one of these 
parameters as t. For the moment, we do not assign to it any physical meaning. 
Of course, the notation anticipates that it will be identified in due course with 
the time variable. 

Let us start by assuming that we are given the divergence p of a field, which 
is a function not only of position but of the parameter t as well, which we now 
make explicit, i.e., p — p(t,r). Let our boundary condition be such that p 
vanishes at spatial infinity. Helmholtz theorem states that there is a field, call 
it D, satisfying 

V-D(f,r)=p(*,r). (30) 
As we have seen, the field D(f, r) is explicitly given by 

D(t, r) = -V r f - P{t ' ri) , d 3 n + V r x Z(t, r), (31) 
J 47r |r — ri | 

with Z an arbitrary field that we are free to put equal to zero, if we want. 
We stress that t plays, in all of this, only the role of a parameter that can be 
appended to the fields, without having any dynamical meaning. The field D(t, r) 
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is required to satisfy only one condition we have put upon it, i.e., V • T>(t, r) = 
p(t, r). The curl of D has been left unspecified, or else set equal to zero. 

Consider now a field j(t, r) depending on the same parameter t as p does. 
Assume next that p(t,r) and j(t, r) satisfy a continuity equation: 

d tP + V-i = Q. (32) 
By using Ea. l|3U[) the continuity equation can be written as 

V-(fitD+j)=0. (33) 

The divergenless vector <9 t D + j can thus be taken as being the curl of a 
field H(t, r). Indeed, according to what we have seen before, the equation 
VxH=j + 9f D can be solved as 

m, r) = V r x / KM^+ftD^) ^ + VrV (t , r) . (34 ) 
7 47r|r-r 1 | 

As long as we do not specify V • H the function V remains undetermined. 
In any case, the Maxwell equations V • D = p and VxH — <3jD = j hold true 
as a consequence of the continuity equation and Helmholtz theorem. However, 
these equations are not enough to determine the dynamics of the fields D and 
H, even though we may ascribe to t the meaning of time. This must be so 
because - to begin with - the continuity equation alone does not entail enough 
information about the dynamics of p and j. But even in case we were provided 
with the complete dynamics of p and j, from a physical point of view it is clear 
that some assumptions must be made concerning the properties of the medium 
(e.g., "space-time") in order to fix the dynamics of the electromagnetic fields 
that will eventually propagate in such a medium. 

Nonetheless, let us pursue a little bit further the mathematical approach 
suggested by Helmholtz theorem. The potentials U and C in terms of which we 
defined the field v(x) read here 



"<*•'> = (35) 

J 47r |r — i"i 

respectively, and we have that D(i, r) = — V r (p(t, r) + V r x Z(t, r) and H(i, r) = 
V r x A(t, r) + V r V{t, r). We obtain then, from Ea. [|35jl . 

For the special case of a point-like charge moving along the curve i"o(£) we 
put p(t, r) = qd(r — ?o(t)) and the above expression reduces to 

D(t,r) = 2 I r ~ r ° ( /!i + v r x Z(t,r). (38) 

47r|r-r (i)| 2 |r-r (t)| 
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According to Ecis. l|37[) or 138|) the field D(t, r) at time t entails an instan- 
taneous Coulomb field produced by a continuous charge distribution p, or else 
by a point-like charge q. Such a result would correspond to an instantaneous 
response of the field to any change suffered by the charge distribution. That 
would be in contradiction with the finite propagation-time needed by any sig- 
nal. Whatever the field Z(i, r) might be, it must contain a similar instantaneous 
contribution that cancels the former one, if we want the present approach to 
bear any physical interpretation. Such an issue has been discussed and cleared, 
in the case of the complete set of Maxwell equations, by showing that both the 
longitudinal and the transverse parts of the electric field contain instantaneous 
contributions, which turn out to cancel each other Note also that by taking 
Z equal to zero in Eg. 14411 we have V x D = in our case, which is not what hap- 
pens when D has to satisfy (together with H) the complete system of Maxwell 
equations. In any event, as we have already stressed, it is necessary to add some 
additional information to the one derived from the continuity equation, in order 
to fix the dynamics of the fields. We do this in the following form. Instead of 
taking the potentials ip and A as given by Eas. (|35|) and l|3tj|) . we assume them 
as additional quantities, out of which we define the fields E and B through 

E(i,r) = -V r p(t,r)-fitA(t,r), (39) 
B(i,r) = V r xA(t,r). (40) 

These fields obey then the homogeneous Maxwell equations identically: 



V-B = (41) 
V x E + d t B = 0. (42) 

Side by side to these two Maxwell equations we write the inhomogeneous 
ones: 



V • D = p (43) 
VxH=j + <9 t D. (44) 

We stress once again that - up to this point - the homogeneous and the 
inhomogeneous equations are independent from one another. We may connect 
them through some constitutive equations, like, e.g., 



D = eE, (45) 
H = fi^B. (46) 

These equations are usually assumed to describe a linear medium of electrical 
permittivity e and magnetic permeability fi. A particular case of such a medium 
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is vacuum, and the system of equations, Eas. (|4"Tl l4*2l 14*51 l4*4*j l. that arises out of 

a connection like the one given by Eas. (|45l I46fl is what we know as the complete 
system of Maxwell equations. 

Without connecting (D, H) with (E, B) through some constitutive equa- 
tions, we have no closed system. The equations that we have written down for 
(D, H), that is Maxwell equations with sources, can also be written down for 
a fluid, for example. We can expect that any conclusion that can be derived in 
the realm of electrodynamics from the equations V • D = p and VxH=j + <9tD 
without coupling them to the source- free Maxwell equations, will have a corre- 
sponding result in the realm of fluid dynamics. This assertion can be illustrated 
by two examples: 1) A fluid having a point-like singularity in its density p 
(one can achieve this approximately, by using an appropriate sink) : one obtains 
in this case a velocity-field obeying a law that is mathematically identical to 
Coulomb's law [Jj. 2) A fluid where a so-called vortex tube appears (tornadoes 
and whirl-pools are associated phenomena) , in which case - after approximating 
the vortex-tube by a line singularity - one obtains a velocity-field through an 
expression which is mathematically identical to the Biot-Savart law 0. 

4 Tensor ial formulation 

The derivation of the inhomogeneous Maxwell equations as a consequence of 
charge conservation is nothing new ^3 ID] ■ ^ follows as a direct application 
of a theorem of de Rahm for differential forms El EOl • According to this 
theorem, given a four- vector j a {x) for which a continuity equation holds, i.e., 
d a j a = 0, there exists an antisymmetric tensor F a @ = —F@ a fulfilling d a F a P = 
j 13 . As we said before, this last equation is nothing but the tensorial form of 
the inhomogeneous Maxwell equations, Eas. (|43|) and l|44|) . Now, the tensor F a @ 
is not always derivable from a vector A a . In order to be derivable from A a in 
the form F afi = d a A - d^A 01 , it must satisfy the equation d a F^ + d^F^ a + 
_ q xhis is the tensorial form of the homogeneous Maxwell equations, 
Eas. (|41l42(l . In other words, given j a and A a , with j a satisfying a continuity 
equation, we may introduce two antisymmetric tensors, F^j and . The first 

one can be determined so as to satisfy d a F^ = j@ , according to de Rahm's 
theorem. The second tensor, defined through F^S = d a A 13 — d@A a , satisfies 

® a F(2) + ® >3 F(2) + d 7 F(2) = identically. In order that these two equations do 
conform a closed system, i.e., the total system of Maxwell equations, we need 
to connect F^ with F^j through some constitutive relation. In the following 
we ellaborate on all this, but without employing the tools of differential forms 
on manifolds, which - in spite of their usefulness - cannot be said yet to be part 
of the lore of physics. 

It is indeed not necessary to rest on de Rham's theorem and the theory of 
differential forms on manifolds, in order to derive the foregoing conclusions in 
tensorial form. One could start with the tensorial form of Helmholtz theorem 
jU El 0] and go-ahead with a similar reasoning as the one we have followed in 
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the preceding sections. We shall however proceed by explicitly constructing a 
tensor fulfilling our requirements. 

Let us thus start by assuming that we are given a vector field j a . We want 
to show that there is an antisymmetric tensor F a @ fulfilling 

d a F^ = /. (47) 

Note first that from Eq. 147( 1 and the antisymmetry of F a P it follows that j 
must satisfy the continuity equation: 

fy/ = °- ( 48 ) 
We now demonstrate the existence of the tensor F a @ by explicit construction. 
To this end, we consider the Green function G(x,x') satisfying 

d li d»G{x,x') = 5{x-x'). (49) 
Given G (x, x') and j a we introduce the potential A^(x) as 

A»(x) = j G(x,x')j»{x')d 4 x', (50) 

and define 

F^(x) = d ft A u (x)-d v A' i (x) (51) 

[&*G (x, x') f(x') - d v G (x, x') f{x')] d A x'. (52) 



Let us now take the four-divergence of the above defined tensor F^ v (x) : 



d ll F» v {x) = J [d^G(x,x')f(x')-d l ,d''G(x,x')j>*(x')}d i x'. (53) 

Because G(x,x') satisfies Ea. (|49|l . the first integral in Eg .1 (53( 1 is equal to 
j v {x). As for the second integral, in order to show that it is zero we do as 
follows. Because G (x, x') satisfies Eg . (|49fl . it must be a function of (x — x'), so 
that <9 M G (x, x') = —d' tl G(x,x'), where gf, = d/dx'^. We use this property and 
integrate by parts the second term in (|53|) ; at the same time we replace the first 
term by j u {x): 



3^{x) = f{x)+d v J[d'^G(x 1 x')f{x'))~G{x,x')d^f{x')]d 4 i&A) 

= f{x) + d v fd' ll (G(x,x')j^x'))d' k x'. (55) 

We may now employ the generalized Gauss theorem in order to show that the 
four- volume integral on the right-hand side of 1(55(1 vanishes. The four- volume 
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has as its boundary a three-dimensional hypersurface S' whose differential ele- 
ment we denote by dS' p . Thus, because j M vanishes at spatial infinity, 

J d' p (G (x, x')j^(x')) d 4 x' = fGix^j^dS^O, (56) 

when we let S' — > oo, and with this result Eq.JSHJ) reduces to 

Now, just as in the three-dimensional case, where the divergence of a field 
did not determine the field uniquely (see Eq. l|28|0 . by subjecting F^ v to the sole 
condition of fulfilling Eq. (|47|l we do not fix completely. Indeed, the tensor 
K 1 *", which is defined below in terms of an arbitrary four- vector B p , fulfills also 
Eq.gJl: 

K nu = piiv _ }_ e fiv P <7 (p p B a — d a B p ) (57) 
= F* v - -e^H^ = F* v - H» v . (58) 

Here, e^ pa is the totally antisymmetric Levi-Civita tensor (in fact, a tensor 
density). The four-divergences of K^ v and F^ v are the same because, due to 
the antisymmetry of e MI/pcr and the symmetry of partial derivatives like dpd p , 
we have 

g^u = 1 e »u P „ ^ dpBa _ d„d a B p ) = 0. (59) 
Hence, we obtain Maxwell equation with sources: 

d a K a P = /, (60) 

together with the identity d p F v \ + d v F\ p + d\F pu = 0, which follows from the 
definition of F» u , as given in Ea.flfl). Introducing the dual F^ u = e^ al3 F a p/2 
of the tensor F a p we can write the former identity in the form 

d p F" u = 0. (61) 

Eas. (|60|l and (|61|l are Maxwell equations in tensorial form. As we said before, 
they constitute a closed system as long as K a @ and F^" become related to each 
other by some constitutive equations. A general, linear algebraic, relationship 
between these tensors has the form 

K af) = l - X af)pa F pai (62) 

where ■y^Pf" 7 is called the constitutive tensor. It has the following symmetry 
properties: ^ a ^P (J = — ^ a P a = „^«/ 5o 'p = i n three-dimensional notation 

the components of K a @ are D and H, whereas those of are E and B. For 
free-space, the only nonzero components of y^P?" have the values eq and l//Jo, 
corresponding to the electrical permittivity £o and magnetic permeability /io of 
the vacuum. 
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The properties of the medium can be specified through an equation like , 
as well as by introducing some other quantities that describe the polarization 
and magnetization of the medium. In three- vector notation these quantities are 
the vectors P and M, respectively. Their relation to (D, H) and (E, B) is given, 
in the simplest case, by 



D 
H 



= £ E + P, 

= — B M. 

Mo 



(63) 
(64) 



In tensor notation, P and M are subsumed into an antisymmetric tensor: 



M 



a/3 



( 


Pi 


P2 


P 3 


\ 


-Pi 





-M s 


M 2 


-Pi 


M 3 





-Mi 




V ~Ps 


-M 2 


Mi 





/ 



(65) 



This choice corresponds to the assignment Ei = for the electric field, and 
Bi = —EijkFjk/l for the magnetic field (Latin indices run from 1 to 3). By relat- 
ing M al3 to K al3 through K a ^ — F a/3 — M a/3 we can rewrite the inhomogeneous 
Maxwell equation (|60|l as 



d a F 



a/3 



■ f + d a M a P. 



(66) 



Written in this form, the inhomogeneous Maxwell equation makes the magnetization- 
polarization tensor M a @ appear as a source of the electromagnetic field F a @ = 
d^A v {x) — d"A^(x). The constitutive equations are given in this case as a 
connection between and F a @ ^31 • At any rate, one has to make some hy- 
pothesis concerning the electromagnetic properties of the medium - be it vacuum 
or any other kind of medium - in order to obtain the closed system of Maxwell 
equations. The simplest assumption is to attribute to the medium the prop- 
erty of reacting locally and instantaneously to the presence of a field. It could 
be, however, that such an assumption describes reality as a first approximation 
only. 

Finally, we want to stress the central role played by the Green function 
G(x, x'). We assumed this function to satisfy Ea. (|49|) . an equation entailing the 
velocity of light. One possible solution of l|49|l is given by the retarded Green 
function 



G(x,x') = 



1 



47r r - 



5{t-t'-\v-v'\/c). 



(67) 



This is the solution of Ea. (|49|l to which we ascribe physical meaning. By using 
it in Ea. H5Ufl we are actually prescribing how the source j^{x') at a space-time 
point x' gives rise to an electromagnetic field A^(x) at a distant point x; a field 
that virtually acts upon a second charge or current density that is located at such 
a distant point. There is therefore a fundamental piece of information concerning 
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the electromagnetic properties of space-time that is already contained in the 
Green function, be it given through the special form of G(x,x'), as in Eq.lJBTJl, 
or through the equation it has to satisfy, e.g., Eq. p^|) . 

5 Summary and Conclusions 

We showed that starting from charge conservation one can arrive at equations 
which are mathematically identical to Maxwell equations with sources. These 
equations are therefore tightly linked to a general statement telling us that 
something is conserved. Consider anything - charge, matter, or whatever - that 
is contained inside an arbitrary volume. Consider also that this thing is in a 
quantity that changes with time. If the change is exclusively due to a flow 
through the volume's boundary, a continuity equation holds true. Then, as a 
consequence of it, a pair of Maxwell-like equations must be fulfilled by some 
auxiliary fields, which take the role ascribed to D and H in Maxwell equations. 

That Maxwell equations with sources follow from charge conservation is a 
mathematical fact that has been known since a couple of years although 
it is usually not mentioned in standard textbooks of electromagnetism. Maxwell 
equations with sources involve the fields D and H, whereas the source- free equa- 
tions involve the fields E and B. It is through some constitutive equations con- 
necting (D, H) with (E, B) that we obtain a closed system, i.e., the complete 
system of Maxwell equations. The constitutive equations express, in some way 
or another, the underlying properties of the medium where the fields act or are 
produced. From this perspective, the Maxwell equations entail besides charge 
conservation some properties of the medium, yet to be unraveled. These prop- 
erties are effectively described, in the simplest case, through the permittivity 
e and the permeability /i of the medium. The first one refers to electrical, the 
second one to magnetic, properties of the medium, be it vacuum or any other 
one. It is just when the equations for (D, H) together with those for (E, B) do 
conform a closed system, that we can derive a wave equation for these fields. 
The velocity of wave propagation is then given by c = (e^) -1 / 2 , the velocity of 
light. This must be in accordance with the assumptions we make when choos- 
ing a physically meaningful Green function. It is remarkable that the velocity 
of light can be decomposed in terms of a product of two independent parame- 
ters. However, the development of physics has led us to see c as a fundamental 
constant of Nature, instead of e and ^. Nevertheless, currently discussed and 
open questions related to accelerated observers, Unruh radiation, self-force on a 
charge, magnetic monopoles and the like, might well require an approach where 
the role of c recedes in favor of quantities like e and fi. Maxwell equations, when 
written in the - by now - most commonly used Gaussian units, do not include 
but the single constant c, hiding so e and \i from our view. These last two con- 
stants might well be key pieces that remain buried under the beauty of a unified 
theory of electromagnetic phenomena, which is the version of electrodynamics 
that we know and use today. A version that should not be regarded as a closed 
chapter in the book of classical physics. 
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